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Jacobi elliptic functions sn(x,m), cn(2;,m) and dn(x,m), with elliptic modulus parameter m = 
/c^ (0 < m < 1) play an important role in describing periodic solutions of many linear and nonlinear 
differential equations of interest in diverse branches of engineering, physics and mathematics. The 
Jacobi elliptic functions are often defined with the help of the elliptic integral 
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Over two centuries ago, John Landen [ffl] studied the quadratic transformation of variables 



t 



This transformation yields the transformed integral 
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These celebrated relations are known as the Landen transformation formulas. They have the special 
property of providing a non-trivial connection between Jacobi elliptic functions with unequal elliptic 
modulus parameters m and m = (1 — \Jl — ■mf' j (1 + \Jl — m)^. 

The purpose of this paper is to give a generalization of all three Landen formulas. To describe our 
approach, let us focus on eq. (^) first. Using the identity dn^(ii, m) = 1 — m sn2(u, m), and changing 
variables to a; = (1 + one can re-write the "dn" Landen formula @ in the alternative form 
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Here, the right hand side contains the sum of two terms with arguments separated by K{m), the 
complete elliptic integral of the first kind Q . Our generalized Landen formulas will have p terms on 
the right hand side. The results are somewhat different depending on whether p is an odd or even 
integer. The generalization of eq. (^) for p odd is given by 



dn(x, m) = Qi ( dn[Qix, m] + dn[aix + 4:K{m) /p, m] + ■ ■ ■ + dn[aix -|- 4{p — l)K{m) /p, m 



(5) 



where 

ai = |^dn[0, m] + dn[4i^ (m) /p, m] + ■ ■ ■ + dn[4(p — l)K{m) /p, m] 
fh = {m-2)al + 2alAi , Ai = ( dn^[0,m]+dn^[4:K{'m)/p,m]-\ \-dn^[A{p-l)K{m)/p,m\] . (6) 



Likewise, our generalization of eq. (y) for p even is given by eq. p^). Similarly, our results for 
generalized "cn" Landen formulas corresponding to eq. (|2|) are eqs. ( p9[ ) (p odd) and (^) {p even) 
and the generalized "sn" Landen formulas corresponding to eq. @) are eqs. (p8|) {p odd) and (p 
even). The richness of the generalized results is noteworthy - most formulas involve sums, but one [eq. 
(|5^] has products; most formulas have all positive signs, but one [eq. (^)] has alternating signs; most 
formulas have non-trivial scalings of the argument u. This large variety of results is a consequence 
of the many different types of periodic solutions for nonlinear equations which we have previously 
obtained |, |]. 

Generalized "dn" Landen Formulas: Given the diversity of the generalized Landen formulas, it 
is necessary to establish them one at a time. To get an idea of our general approach, let us first focus 
on the proof of eq. (^. Consider the periodic solutions of the static sine-Gordon field theory in one 
space and one time dimension, that is, the solutions of 

4>xx = sin . (7) 

Note that the time dependent solutions are easily obtained from here by Lorentz boosting. One of the 
simplest periodic solutions of the field eq. <^ is given by 

sin(</)(x)/2) =dn(x,m) . (8) 

It was shown in refs. [^, |5| that a kind of superposition principle works even for such nonlinear 
equations because of several highly nontrivial, new identities satisfied by Jacobi elliptic functions [^. 
In particular, it was shown in ref. [^ that for any odd integer p, one has static periodic solutions of 
eq. (^) given by 

p 

sin((/)(x)/2) = ai ^ , p odd , (9) 



i=l 



where 



with ai given by eq. 
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The question one would like to address here is if solution (^) is completely new, or whether it can 
be re-expressed in terms of simpler solutions like (^), but where m and rh need not be the same. To 
that end, we note that on integrating the field eq. (|^) once, we obtain 



bl = C -2cos(j) . (11) 



On further integration, this yields 



x + xo, (12) 



"v/C — 2 cos 1 

where xq is a constant of integration which we put equal to zero without loss of generality. On 
substituting 

sin((/./2) = ^ , (13) 



equation (12) takes the form 



X . (14) 



Now the important point to note is that if we perform the integral for different values of C then 
we will get all the solutions of eq. (0). Further, if two solutions have the same value of C, then they 
must be the same. As far as the integral ([l^ ) is concerned, it is easily checked that the three simplest 
solutions covering the entire allowed range of C are 

tp = sech X , C = 2 , (15) 

ip = dn(x, m) , C = Am — 2 , (16) 

V' = cn(^,m), C= — -2, (17) 

where < m < 1. Note that the constant C has been computed here by using eq. (0), which m 
terms of ^(x) takes the form 

C = 2 - 4^2 + . (18) 

1 — 

Thus, whereas for the solution ([l^), C lies in the range —2 < C < 2, for the solution ([T7|) , C lies 
between 2 and co. Note that for C < — 2, there is no real solution to eq. (jl^). 

Now the strategy is clear. We will take the solution (^) and compute C for it and thereby try to 
relate it to one of the basic solutions as given by eqs. (|l^) to ([T7|). One simple way of obtaining the 
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constant C from eq. (|18|) is to evaluate it at a convenient value of x, say x = 0. In this way, we find 
that for the solution (^), C is given by 

C = -2 + 4(m - 2)al + 8a?^i , (19) 

where ai and Ai are as given by eq. (P). Now, as m ^ 0, ai = dn(x,m = 0) = 1 and hence 
C = —2. On the other hand, as m ^ 1, K{m = 1) = oo, dn(x,m = 1) = sech x and hence ai = 1 
so that C = 2. Thus for solution (|9|), as m varies in the range < m < 1, the value of C varies in 
the range — 2 < C < 2. Hence it is clear that the solutions (P) and (|l^) must be same. On equating 
the two C values as given by eqs. (|l^) and ([l^), we find that the two solutions are identical provided 
m and m are related by eq. (|6|) and hence the appropriate Landen transformation valid for any odd 
integer p is given by eq. (^) with m and fh being related by eq. (P). 

What about the case of even p? We have checked Q that for this case, an exact solution is 

p 

sin(</)(x)/2) = 02 ^ , p even (20) 
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where 
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with 02 being given by 

a2 = (^dn[0, m] + dn[2K{m) /p,m]-\ h dn[2{p - l)K{m) /p, m]J . (22) 

Proceeding exactly as before, we find that C is again in the range —2 < C < 2, and hence comparing 
with solution ( [l^ ) yields the Landen transformation for even p: 

dn(x, m) = 02 ^dn[o2X, m] + dn[o2X + 2K{m) /p,m] + ■ ■ ■ + dn[o2X + 2{p — l)K{m) /p, m]j , (23) 

where 02 is given by eq. ( p^ and rh is given by 

m = (m- 2)03 + 20^^2 , A2 = ( dn^[0,m]+dn^[2/i:(m)/p,m]H hdn^[2(p- l)K(m]/p, m) ) . (24) 



Note that when p = 2^ one recovers the Landen formula (^), since dn [K{m),m] = k', and rh simplifies 

to {i-k')y{i + k')\ 

Generalized "cn" Landen Formulas: As shown in ref. j^, another periodic solution of the static 
sine-Gordon eq. @j is 

p 

sin((/)(x)/2) = os^Ci , p odd , (25) 



i=l 



where 



with Q3 being given by 



a^x 4:{i — l)K{m) 

Cj = cn(^= H ,m) , (26) 

^/m p 



03 = ^cn[0, m] + cn[4K{m) /p, m] + • • • + cn[4(p — l)K{m) /p, m]^ . (27) 
Using eq. ([l8| ) we can now compute the corresponding value of C. We obtain 

C = -2+-^ ^+8ai^3 , A3 = icn^[0,m]+cn^[AK{m)/p,rn] + ---+cn^[A{p-l)K{rn]/p,rn)j . 

(28) 

It is easily checked that since < m < 1, C varies from 2 to 00, and hence the solutions ( [T7| ) and (|25|) 
must be identical. On equating the two values of C as given by eqs. (|l^) and (^8[), we then find that 
for odd p, the Landen transformation is given by 

cn(x, m) = as ^cn [/?x, m] + cn + AK{m) /p, m] + • • • + cn [f3x + 4:{p — l)K{m) /p, m] ^ , (29) 

where 03 is given by eq. (p7|), and /3 and m are given by 



' ~ (1 - 2m)a| + 2ma|^3 ' ^ 



with ^3 being given by eq. (28). 

What happens if p is an even integer? As shown in |^] , in that case another periodic solution of 
the static sine-Gordon eq. (^) is 

p 

sin((/)(x)/2) = 04 ^ [di - dj+i] , p even , (31) 

i odd 



where di is as defined by eq. (|2lD while 04 is given by 

04 = ^dn[0, m] — dn[2K{m) /p,m] + ■ ■ ■ — dn[2{p — l)K{m) /p, m]^ . (32) 

Using eq. (|^) the value of C for this solution is easily computed and we find that 2 < C < cxd. On 
comparing with solution (|l^ we find that in this case the Landen formula is 

cn(x, rh) = 04 ^nla^V^x, m] — dn[a4\/?fix + 2K{m) /p, m] + • • • — dn[a4,Vmx + 2{p — l)K{m) /p, m]J , 

(33) 

where 04 is given by eq. (p^ while rh is given by 



171 = 3—, A4 = (dn^[0,m]-dn^[2K(m)/p,m]H dn^[2{p-l)K{m)/p,m]) . (34) 



(m - 2)q| + 2a|^4 
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Generalized "sn" Landen Formulas: Here, we start from the sine-Gordon field equation (c = 1) 

(t>xx - 4>tt = sin (/) , (35) 

and look for time-dependent solutions with velocity v > 1 (which are called optical soliton solutions 
in the context of condensed matter physics). In terms of the variable 

X — vt , , 

V = -=^= , (36) 
Vf — 1 



the field eq. (35) takes the form 

(pr,ri = - sin . (37) 
On integrating this equation once, we obtain 

(P'^ = C + 2cos(p . (38) 

On integrating further, we get 

/ /n ft A. = ^ + ^0 , (39) 
J \JC -I- 2 cos 

where r/o is a constant of integration which we put equal to zero without loss of generality. Substituting 
sin((/)/2) = V', yields 

T] . (40) 



If we now perform the integral for different values of C, then we get all the solutions. It is easily 
checked that the three simplest solutions of eq. (^) covering the entire allowed range of C are 

'0 = tanhr/, C = 2, (41) 

V' = \/?fisn (?7, m) , C = 4m — 2 , (42) 
V' = sn(^,7fi) , C=l-2, (43) 



where < m < 1. Note that the constant C has been computed here using eq. (38), which in terms 
of 'il) takes the form 

C7 = -2 + 4V^2 ^ . (44) 

\ — '^)^ 

Thus, for solution (p2|), C is in the range —2 < C < 2, whereas for solution (|4^), C lies between 2 and 
CO. Note that for C < —2, there is no real solution to eq. (^0|). 
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Using appropriate linear superposition, it was shown in ref. ||5| that for odd p one of the solutions 
of eq. (37) is given by 

p 

sn\{(f){vi) /2) = v^ai ^ Sj , p odd , (45) 



where 

with ai being given by eq. 
We find 



Si = sn 



A(i - l)K(m) 

aiT] H ,m 

P 



(46) 



Using eq. ([44D , we can now compute the corresponding value of C. 



C = -2 + 



a 



2 ' 



(47) 



where ai,a3 are given by eqs. (0) and (p7|) respectively. It is easily checked that since < m < 1, C 



has values between -2 and 2. Hence the solutions (42) and (|45| ) must be identical. On equating the 
two values of C as given by eqs. ( ^2] ) and (^7|), we then find that the "sn" Landen transformation 
formula for odd p is given by 



sn(x, m) = Q3 ( sn[aix, m\ + sn[aix + 4i^ (m) /p, m] + • • • + sn[aix + 4(p — l)K{m) /p, 



m, 



with m given by 



m = m- 



a 



2 • 



(48) 



(49) 



Finally, we turn to the "sn" Landen transformation formula for the case when p is an even integer. 
One can show ||] that in this case, a solution to eq. ( |37| ) is given by 

(50) 



sin(0(?7)/2) = mF^'^a2A^Ii^^^^Si , p even 



where 



Si = sn 



2(i - l)K{m) ■ 
a2'q -I , m 



P 



with a2 being given by eq. (p2|) and defined by 



= sn[2i^ (m) /p, m] sn[AK{m) /p, m] • • • sn[2(p — l)K{m) /p, m] 



Using eq. (|44|), we can now compute the corresponding value of C. We obtain 



(51) 



(52) 



C = -2 + AmPa^Ai . 



(53) 
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It is easily checked that since < m < 1, the value of C varies between -2 and 2 and hence the 
solutions (42) and ( |50| ) must be identical. On equating the two values of C as given by eqs. (42) and 



(p^, we then find that for even p, the Landen transformation formula is 

A5a2sn(x, in) = SD.[a2X^ m] sii[a2X + 2K{m) /p, m] • • • sn[a2X + 2{p — l)K{m) /p, m] , (54) 
with m given by 

rh = mFa^A^ . (55) 

It is amusing to notice that as m ^ 0, 

p 

A^{p, m = 0) = sin(7r/p) sin(27r/p) • • • sin((p - l)7r/p) = -—^ . (56) 

Transformed Modulus Parameters: At this point, we have generalized all three of the celebrated 
two hundred year old p = 2 Landen formulas [eqs. (Q), (|2[), (^] to arbitrary values of the gener- 
alization being different depending on whether p is an even or odd integer. One might wonder that 
whereas the relation between rh and m is the same for all three (p = 2) Landen identities, it seems to 
be different for higher values of p. However, quite remarkably, we have established analytically that 
both for p = 3 and for p = 4, the relation between m and m is in fact the same for all three generalized 
identities. For example, we find that for p = 3, the relation between rh and m is the same for all 
three Landen transformations as given by eqs. (P), ( |30|) and (|49|). All expressions can be algebraically 
simplified and written in the common form 

m = m-. -777 , 57) 

(l + g)2(l + 2g)2 ' ^ ^ 

where q = dia{2K{m)/3,m). Note that while deriving this result, use has been made of the fact 
that cn(4:K(m)/3,m) = "j^^ and that q satisfies the identity + 2q^ — 2(1 — m)q — (1 — m) = 0. 
Similarly, using the relations dn{K{m)/2,m) = dn{3K{m)/2,m) = t and dn(i^(m),m) = t^, where 
t = (1 — m)^/^, it is easily proved that for p = A the relation between rh and m is the same for all 
three Landen transformations as given by eqs. (p^), (p^) and (|55D, i.e. m = (1 — t)^/(l + 1)^ . 

What about the results for higher values of pi Although the analytic proof seems rather compli- 
cated, we have nevertheless checked numerically using the mathematical software package Maple, that 
the relation between fh and m is the same for all three Landen formulas with the same value of p. The 
numerical results for m as a function of m for various values of p ranging from 2 to 7 are shown in 
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Table 1. Note that for any fixed value of p, as m increases from to 1, m also increases monotonically 
from to 1. Also, for any given fixed value of m, m decreases monotonically as p increases. 

In all the Landen formulas (|^), (|2|) and the modulus parameter rh is less than the modulus 
parameter m, and this is often called an ascending Landen transformation j^]. Note that the gener- 
alizations which we have established in this paper in terms of shifts involving the period K{m) on 
the real axis maintain the relationship m < m. We are currently working on the opposite case of de- 
scending Landen transformations and their connection with additional generalized formulas in terms 
of shifts involving the period K'{m) on the imaginary axis. 

The results of this paper clarify the relationship between the well known periodic solutions of the 
various nonlinear equations and those obtained by us using the idea of judicious linear superposition 
^, |5|. Further it provides a deep connection between the highly nonlinear Landen transformation 
formulas involving changes of the modulus parameter and certain linear superpositions of an arbitrary 
number of Jacobi elliptic functions. 

We are very grateful to the U.S. Department of Energy for providing partial support of this research 
under grant DOE FG02-84ER40173. 
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Table 1: A table of the modified modulus parameter fh in the generalized Landen transformation 
formulas as a function of the modulus parameter m and the number of terms p in the formula. Note 
that for odd integers p, the values of m are obtained from eqs. (|^, ( |30[ ) or (|49|), and as mentioned in 
the text, they are all the same. Similarly, for even integers p, the values of m are obtained from eqs. 
(p^), (p^) or (|55D, and they are also all the same. 
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rh{p = 2) 


■m{p = 3) 




rh{p = 4) 




m{p = 5) 




m{p = 6) 




■m{p = 7) 



































0.25 


.5155 X 10-2 


.9288 X 10" 


-4 


.1669 X 10" 


-5 


.3000 X 10" 


-7 


.5392 X 10" 


-9 


.9693 X 10" 


11 


0.5 


.2944 X 10"^ 


.1290 X 10" 


-2 


.5580 X 10" 


-4 


.2411 X 10" 


-5 


.1042 X 10" 


-6 


.4503 X 10" 


-8 


0.75 


.1111 


.1005 X 10" 


-1 


.8666 X 10" 


3 


.7438 X 10" 


-4 


.6381 X 10" 


-5 


.5475 X 10" 


-6 


0.9 


.2699 


.4311 X 10" 


-1 


.6158 X 10" 


-2 


.8655 X 10" 


3 


.1213 X 10" 


3 


.1701 X 10" 


-4 


0.99 


.6694 


.2506 




.7283 X 10" 


-1 


.1963 X 10" 


-1 


.5185 X 10" 


-2 


.1362 X 10" 


-2 


0.999 


.8811 


.5292 




.2374 




.9312 X 10" 


-1 


.3464 X 10" 


1 


.1264 X 10" 


-1 


0.9999 


.9608 


.7446 




.4481 




.2293 




.1080 




.4891 X 10" 


-1 


0.99999 


.9874 


.8721 




.6374 




.3973 




.2239 




.1193 
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1 


1 




1 




1 




1 




1 
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